Enhancement of the superconducting transition temperature in Nb/Permalloy 
bilayers by controlhng the domain state of the ferromagnet. 
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In (S/F) hybrids the suppression of superconductivity by the exchange field hex of the ferromagnet 
can be partially lifted when different directions of hex are sampled simultaneously by the Cooper 
pair. In F/S/F trilayer geometries where the magnetization directions of the two F-layers can be 
controlled separately, this leads to the so-called spin switch. Here we show that domain walls in a 
single F-layer yield a similar effect. We study the transport properties of Nio.80Feo.20/Nb bilayers 
structured in strips of different sizes. For large samples a clear enhancement of superconductivity 
takes place in the resistive transition, in the very narrow field range (order of 0.5 mT) where 
the magnetization of the Py layer switches and many domains are present. This effect is absent 
in microstructured samples. Comparison of domain wall width 5^ to the temperature dependent 
superconductor coherence length £,s{T) shows that 5^ ~ ^s{T), which means that the Cooper pairs 
sample a large range of different magnetization directions. 

PACS numbers: 74.45.+c,74.78.-w,85.25.Hv 



Proximity effects between a superconductor (S) and a 
ferromagnet (F) are the focus of much current research, 
basically because of the possibilities for several distinct 
and unusual phenomena. One is due to the fact the ex- 
change field hex in the ferromagnet gives rise to an oscil- 
latory damped amplitude of the pairing function. In an 
S/F/S geometry, this oscillation allows coupling of the 
superconducting banks with a phase change of tt rather 
than P, Q ■ This can be witnessed in a nonmonotonic 
variation of the superconducting transition temperature 
Tc as function of the F-layer thickness dp- 0,01; or, espe- 
cially with a weak ferromagnet, in the Josephson current 
of an S/F/S junction 5, 6] or a superconducting quan- 
tum interference device (3 • Other phenomena are linked 
to a situation in which the Cooper pair can sample dif- 
ferent directions of hex within its coherent volume. The 
best known example is an F/S/F geometry with a thin 
S-layer, in which the magnetization of one F-layer can 
be rotated with respect to the other. The suppression of 
the order parameter in S will then be larger when both 
magnetizations are parallel (P) and smaller when they 
are antiparallel (AP). This so-called spin switch has been 
discussed theoretically 0, and a first experimental re- 
alization was recently reported by Gu et al. They 
used a weak ferromagnet (CuNi) and reported a small 
but measurable increase of with respect to Tf (the 
transition temperatures of the P and AP configurations, 
respectively) of about 5 mK. However, the physics of the 
problem is more general. A weak ferromagnet is not an 
a priori condition for the effect, and it might even be ar- 
gued (although this has not been emphasized in theoret- 
ical treatments) that ferromagnets with large exchange 
fields are preferable. Also, the trilayer configuration is 



not the only one which can invoke different directions 
for hex '■ any domain wall in the ferromagnet offers dif- 
ferent directions intrinsically, both in the wall and on 
either side. In claiming a coupling effect between two 
F-layers, it might even be necessary to check the absence 
of in-plane domain wall effects. Moreover, other mech- 
anisms may be at play; inhomogeneous exchange fields 
are predicted to induce enhanced superconductivity by 
spin-triplet excitations [ill H^. 

In this Letter, we investigate S/F bilayers and show that 
the domain state of a magnet with strong spin polar- 
ization (Ni8oFe2o, Permalloy (Py)) gives rise to an even 
slightly larger increase in Tc of the superconductor, in 
our case Nb, than mentioned above. We give a qualitative 
discussion of the possible mechanism for the effect, which 
is basically due to a lowered average exchange field seen 
by the Cooper pair and might either be called an in-plane 
spin switch or domain wall induced Tc-enhancement. 
Samples of Nb/Py were prepared by sputter deposition 
in an ultrahigh vacuum system. They were structured 
in simple bars of either 0.5 mm x 4 mm ('large' sample) 
or 1.5 /zm x 20 iim ('small' sample). Contacts were not 
included in the geometry in order to minimize problems 
with stray fields from contact pads or arms. Instead, Au 
contacts for measuring in 4-point geometry were added 
by sputter deposition. The choice of Py as ferromagnet 
is dictated by the wish for large hex (the spin polariza- 
tion is 45 % [13)1 but equally by the need of well-defined 
magnetization switching at low fields. Since an easy axis 
for magnetization is induced by the residual magnetic 
fields in the sputtering machine, care was taken to align 
the long axis of the bars with the easy axis of magnetiza- 
tion Bg. Magnetic fields were applied in the plane of the 
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FIG. 1: Resistance R versus temperature T of the 'large' 
sample s/Py(20)/Nb(21) in zero field (open symbols) and in a 
field of 4.2 mT (filled symbols). The inset shows the difference 
between the two data sets. 



sample, along the bars and therefore along e^. Different 
layer thicknesses were used for both Nb and Py, which 
yielded similar results. We will concentrate on samples 
with both the Nb and the Py thickness around 20 nm. 
The zero-field resistance of a 'large' sample of 21 
nm Nb on 20 nm Py on a Si substrate (denoted 
s/Py(20)/Nb(21)) as function of temperature is shown 
in Fig. n The transition temperature is around 5.7 K, 
depressed from the pure Nb value by the proximity of the 
F-layer. The width is about 100 mK. After stabilizing the 
temperature at around 5.65 K, in the transition region, 
the magnetic field noHa is swept from 65 mT down to 
-65 mT and back. The behavior of R{Ha) is shown in 
Fig. Efc). From the positive field side at 15 mT, R ini- 
tially goes down, shows a small increase around 5 mT, 
but then goes down again, followed by a steep dip in a 
very small field regime around Hdip = -4.2 mT. Revers- 
ing the sweep from the negative field side, the behavior 
is symmetric, with a resistance dip now at -1-4.2 mT. 
Also shown in Fig. El is the magnetization of the sam- 
ple at 8 K, normalized by the saturation magnetization 
Ms- The loop width is about 8 mT and the switch is 
quite sharp, although some rounding can be seen close 
to the coercive field which is due to the misalignment of 
the easy axis by a few degrees. It is clear that the deep 
dips in R{H) occur precisely at the switching field of the 
magnetic layer, which indicates that the domain state of 
the ferromagnet is involved, which is most prominent in 
the steep part of the magnetization reversal. Since the 
dip is very pronounced, we can also measure i?(T, Hdip) 
and compare this with i?(T, 0) in Fig. We find that 
R{T, Hdip) lies consistently below _R(T, 0), with a max- 
imum difference of about 10 mK. Several samples were 
measured which all show the same effect, and one data 
set for a sample s/Py(20)/Nb(19) (reversed Nb and Py) 
is shown in Fig. 13 Again the deep dips in R{Ha) (inset 
in Fig. 121) are clearly present, but there are differences in 
the details. For instance, Tc is somewhat lower, around 
4.85 K, presumably because Nb layer is slightly thinner. 
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FIG. 2: Left-hand scale : Resistance R versus applied field 
HoHa of the 'large' sample s/Py(20)/Nb(21). Filled symbols 
are in positive field (forward) direction, open symbols in back- 
ward direction. Right-hand scale : magnetization M (dotted 
lines) normalized on the saturation magnetization Ms versus 
fj,oHa, measured at 8 K. In both cases, Ha\\ Be (the easy axis 
of magnetization) . 



Also, the transition is sharper, with a width of 50 mK. 
This necessitates a temperature stability during the field 
sweep of better than 1 mK. The magnetization loop (not 
shown) is also wider, about 14 mT, and the rise in R{Ha) 
before the dip is reached is higher as well. We believe 
both facts may be due to the fact that the rotation of the 
magnetization in the covered Py-layer proceeds slightly 
differently than in the free but probably oxidized Py- 
layer. Again, R{T,Hdip) lies consistently below R{T,0), 
with a smooth behavior of the difference. 
The basic explanation we want to offer for our observa- 
tions is that domain walls are formed in the Py-layer dur- 
ing the switching of the magnetization, and that these do- 
main walls lead to the enhanced superconductivity. In or- 
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FIG. 3: Resistance R vs. temperature T for a 'large' sample 
s/Nb(19)/Py(20) in zero field (open symbols) and in a field 
of 6.3 mT (filled symbols). Dotted lines are guides to the eye, 
the drawn line is the difference between the two data sets. 
Inset : R as function of applied field fioHa- Filled symbols 
are in forward direction, open symbols in backward direction. 
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FIG. 4: Resistance R vs. magnetic field fioHa (a) for a 
'small' sample s/Nb(17)/Py(20); (b) same for a large sam- 
ple s/Nb(20)/A10x(8)/Py(20). Filled symbols are in forward 
direction, open symbols in backward direction. 
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der to argue this better, we performed measurements on 
two other types of samples. Fig.^a,b) show R{Ha) data 
(in the transition) on a 'small' (1.5 /^m x 20 /xm) sample 
s/Nb(17)/Py(20) where no effect is found. The other is 
a 'large' sample s/Nb(20)/A10x(8)/Py(20), where a thin 
Al layer was deposited and oxidized before the Py layer 
was grown. No dips are found, but instead a small resis- 
tance increase is seen in the field region of the Py loop, 
probably due to the effect of stray fields from the mag- 
netic layer on the superconducting layer. This shows that 
a proximity coupling is necessary for the dips to occur. 
Before discussing the effects of domain walls on the su- 
perconductivity in more detail, it is necessary to address 
the characteristics of the walls. For the film thicknesses 
dp use (around 20 nm), they are beheved to be of 
Neel-type : the rotation of the magnetization occurs in 
the plane of the sample rather than out-of-plane (Bloch- 
wall) , with the transition between the two at thicknesses 
around 40 nm 0, 0| . We note that the magnetic flux 
coming out of a Bloch wall would suppress the supercon- 
ductivity rather than do the opposite. The other impor- 
tant parameter is the width Sw of the walls, which in Py 
is large due to the small magnetic anisotropy. For Bloch 
walls, 5w is of the order of 0.5 /im, for Neel walls and the 
case that dp << Sw it is of similar magnitude. A detailed 
study by scanning electron microscopy with polarization 
analysis on thick Py films with surface Neel walls yielded 
0.25 /im for the half-width of the walls [l^. The fact 
that Sw is much larger than the low-temperature coher- 
ence length ^5 of the superconductor (for our Nb, ~ 
12 nm) will come back in the discussion. The absence of 
an effect for the small samples s/Nb(17)/Py(20) can be 
explained by the fact that no stable domains were formed 
during switching of the magnetization direction. To con- 
firm this assumption, we simulated the switching behav- 
ior for a range of Py structures with thickness 20 nm and 



FIG. 5: The switching behavior of small Py structures. The 
top picture represents magnetization vs. external magnetic 
field for lOfim long Py structures with aspect ratio a=l (open 
squares, narrow loop) and a=10 (closed squares, wide loop). 
The bottom picture shows the correspondent domain configu- 
ration near the coercive field for both structures. For samples 
with the aspect ratio above 2.5 we found no stable multi- 
domain state during magnetization switching. 

different size and aspect ratio, using the OOMMFF code 
[l7| |. Values for the saturation magnetization Ms and 
the magnetocrystalline anisotropy Ki were determined 
from the magnetization measurements and taken to be 
Ms = 860x10^ [A/m] and Ki = 500 [J/m^]. The easy 
axis direction was taken 6.5° away from the long axis of 
the structure in order to take the slight sample misalign- 
ment into account . Some pertinent results are shown 
in Fig. [S| In a structure with length £ and width w of 
both 10 /xm, the magnetization loop is small and near the 
coercive field He (around 2 mT) stable domain configu- 
rations are present. Upon decreasing w and increasing 
the aspect ratio a = £/w, stable configurations no longer 
occur above a = 2.5. For the 1 fim x 10 /xm structure, 
without a stable domain configuration, has increased 
to 10 mT. Note that this is still below the value expected 
for a uniform rotation of the magnetization, which is 
of order Ms/2 according to the Stoner-Wohlfarth model 
|l4| , since domains are formed during the switching. The 
point is that they are not stable. Simulations on struc- 
tures with larger length (above 50 ^m) showed stable 
multi-domain configurations near He even for structures 
with a > 10. It is therefore the combination of small 
dimensions and large aspect ratio in our 'small' samples 
which precludes stable domains. 

The mechanism we believe to be responsible for the 
enhancement of the superconductivity in the 'large' 
samples is basically that the pair breaking experienced 
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by a Cooper pair is smaller when it samples different di- 
rections of the exchange field. In a sense, this is the same 
mechanism as responsible for the F/S/F spin switch, but 
the simultaneous sampling of two F-layers yields a dif- 
ferent type of averaging. In order to estimate the effect 
of the domain wall, we first consider the one-dimensional 
case of step-like variation of hex which induces changes 
in the order parameter over a distance (^p (the coherence 
length in the F-layer) leading to a lowered pair breaking 
parameter in the vicinity of the domain wall [l^ and 
superconductivity which is enhanced with respect to 
the depression of the uniform F-layer. The situation 
resembles Tc-enhancement by twin planes |l9t Eoj ] and 
can be treated accordingly. An estimate for the order 
of magnitude of the enhanced critical temperature Tew 
can be made as follows. The variation in pair breaking 
occurring over a distance induces a superconducting 
order parameter over a distance S,s{Tcw) around the wall. 
In that case, the effective change in pair-breaking will be 
S,F /S,s{Tcw), and the Tc-enhancement is correspondingly 
a = {Tcw--Tcf)/Tcf ~ £,F/^siTcw), with Tcf the critical 
temperature due to the homogeneous suppression of the 
F-layer. Taking into account that £,s{Tcw) ~ Cs(0)/\/aj 
we obtain a w (Cf/Cs(0))2. With ^5(0) w 12 nm and 
£,F for the strongly magnetic Py « 1 nm, the effect 
is in the range of 1 %. The assumption of a step-like 
change in hex is not correct, given the large value for 
Sw, but it gives a feeling for the orders of magnitude. 
For larger (5^,, the effect will increase as long as 
< Cs(0), but for still larger it has to decrease to zero 
again : the Cooper pair cannot sense the variation in 
hex any longer. In this limit, the relative increase of Tc 
is of order (Cs(0)/^u))^ ^21], which is again about 1 %. 
Another way of arguing that the experiment is sensitive 
to domain wall formation comes from considering 
the temperature dependent (Ginzburg-Landau) £,siT) 
directly. Since we are measuring very close to is 
actually much larger than the low-temperature value. 
From a fluctuation analysis of R{T) we estimate Tc to be 
close to the top of the transition. In Fig. O it would be 
at 4.83 K, where the 0-field and in-field curves start to 
separate. With a typical transition width of 30 niK, the 
relative temperature at the zero of resistance tr = T/Tc 
is 6 X 10"'^, which makes ^s{tr) = S,s{0)/V^ — tr about 
0.15 fim. In the transition, therefore, the condition is 
^s{T) ~ dw, and the Cooper pair samples a considerable 
part of the rotation of the magnetization. This implies 



that in the case of Py the in-plane switch will only 
be visible close to Tc, since at lower temperature the 
magnetization is homogeneous on the scale of ^5. It is 
of interest to note that enhancement of critical currents 
below Tc has been reported for the case of Nb/Co, where 
the domain walls are considerably smaller |23| 
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